Non-Gaussian statistics and extreme waves in a nonlinear optical cavity 
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A unidirectional optical oscillator is built by using a liquid crystal light-valve that couples a pump 
beam with the modes of a nearly spherical cavity. For sufficiently high pump intensity, the cavity 
field presents a complex spatio-temporal dynamics, accompanied by the emission of extreme waves 
and large deviations from the Gaussian statistics. We identify a mechanism of spatial symmetry 
breaking, due to a hypercycle-type amplification through the nonlocal coupling of the cavity field. 

PACS numbers: 42.50.Gy 42.70.Df, 42.65.Hw 77.22. Gm 



Extreme waves are anomalously large amplitude phe- 
nomena developing suddenly out of normal waves, living 
for a short time and appearing erratically with a small 
probability. These rare and extreme events have been 
observed since long time on the ocean surfaces and, 
in this context, have been called freak or rogue waves. 
Recently, rogue waves have been reported in an opti- 
cal experiment Q and in acoustic turbulence in He II 
Q. Associated with extreme waves are L-shaped statis- 
tics, with a probability of large peak occurrence much 
larger than predicted by Gaussian statistics. Different 
mechanisms to explain the origin of extreme waves have 
been proposed, including nonlinear focusing via modu- 
lational instability (MI) [4J and focusing of line currents 
In optical fibers, numerical simulations of the non- 
linear Schrodinger equation (NLSE) [g[ have established 
a direct analogy between optical and water rogue waves. 
In a spatially extended system, the formation of large 
amplitude and localized pulses, so-called optical needles, 
has been evidenced by numerical simulations for trans- 
parent media with saturating self-focusing nonlinearity 
7] . Similar space-time phenomena, such as collapsing fil- 
aments, are also predicted in optical wave turbulence [8J. 
Recently, the spatio-temporal dynamics of Mi-induced 
bright optical spots was observed in Ref. @ and an alge- 
braic power spectrum tail was reported due to momen- 
tum cascade [10|. Nevertheless, no experimental evidence 
has been given up to now of extreme waves in a spatially 
extended optical system. 

Here, we report, what is, at our knowledge, the first ex- 
perimental evidence of extreme waves and non-Gaussian 
statistics in a 2D spatially extended optical system. The 
experiment consists of a nonlinear optical cavity, formed 
by a unidirectional ring oscillator with a liquid crystal 
light-valve (LCLV) as the gain medium. While for low 
pump the amplitude follows a Gaussian statistics, for 
sufficiently high pump we observe large deviations from 
Gaussianity, accompanied by the emission of extreme 
waves that appear on the transverse profile of the op- 
tical beam as genuine spatiotemporal phenomena, devel- 



oping erratically in time and in space. The observations 
are confirmed by numerical simulation of the full model 
equations. Moreover, by introducing a mean-field sim- 
plified model, we show that extreme waves in the cavity 
are generated by a novel mechanism, which is based on a 
hypercycle-type amplification [ill ] occurring via nonlocal 
coupling of different spatial regions. 
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FIG. 1: (color online), a) Instantaneous experimental pro- 
file of the transverse intensity distribution; b) a 1-D profile 
showing an extreme event; / is measured in gray values, g.v. 

The experimental setup is essentially the one de- 
scribed in [12j. The ring cavity is formed by three high- 
reflectivity dielectric mirrors and a lens of / = 70 cm fo- 
cal length. The total cavity length is L = 273.3 cm and 
the lens is positioned at a distance L\ = 88.1 cm from 
the entrance plane of the LCLV. The coordinate system 
is taken such that z is along the cavity axis and x, y 
are on the transverse plane. A LCLV supplies the gain 
through a two-wave mixing process [l3| that couples the 
pump beam with the cavity modes. The LCLV is made 
of a nematic liquid crystal layer, thickness d = 14 /im, 
inserted in between a glass wall and a thin slice 20x30x1 
mm 3 of the photoconductive -B12SW20 (BSO) crystal. A 
voltage Vq is applied by means of transparent electrodes. 
The working point is fixed at Vb — 20.3 V, frequency 75 
Hz. The LCLV is pumped by an enlarged and collimated 
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(10 mm diameter) beam from a solid state diode pumped 
laser (A = 532 nm), linearly polarized in the same direc- 
tion of the liquid crystal nematic director. The pump and 
the cavity field, E p and, respectively, E c , are polarized 
in the same direction and have a frequency difference of 
a few Hz, the detuning being selected by the voltage Vo 
applied to the LCLV Q< 
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FIG. 2: (color online) Experimental PDF of the cavity field 
intensity; the pump is varied from Ip/Ith = 1.8 (red), 4.0 
(blue), 4.2 (green) and 6.4 (purple). The red distribution 
is practically exponential. The black line is the stretched 
exponential function that fits the purple distribution, with 
ci = 29.5. The inset shows the L-shape of the PDF in the 
linear scale. 

A small fraction (4%) of the cavity field is extracted 
by a beam sampler and sent to a CCD camera (1024x768 
pixels and 12 bits depth). The detection optical path is 
set so that the camera records the same field distribution 
as in the entrance plane of the LCLV. The Fresnel num- 
ber is fixed to F = 130 and can be carefully controlled 
through a spatial filter introduced near the lens. 

For relatively low pump intensity, I p = \E p \ 2 = 
2.0 mW/cm 2 , the cavity field shows a complex spatio- 
temporal dynamics, with the formation of many uncorre- 
cted domains. At larger pump, extreme events appear as 
large amplitude peaks, standing over a speckle-like back- 
ground. An instantaneous experimental profile of the 
transverse intensity distribution for I p = 4.2 mW/cm 2 , 
is shown in FigJT] In the inset the 1-D profile shows a 
large event. Note the difference of the peak amplitude 
with respect to the background values of the intensity. 
The locations of the large peaks change spontaneously in 
time. The typical time-scale of the dynamical evolution 
is 100 ms, which is ruled by the response time of the 
liquid crystals (l5j |. 

The probability density functions (PDF) of the light 
intensity are determined experimentally by acquiring a 
large set of images (about one thousand) and then per- 
forming the histograms of the intensity values on the 
whole image stack. In FigJ5]the PDF of the cavity field 



intensity, / = \E C \ 2 : are displayed for different values 
of the pump, I p /I th = 1.8, 4.0, 4.2 and 6.4, I th = 1.2 
mW/cm 2 being the threshold for optical oscillations. The 
dark probability P(I = 0) is the same for all the pump in- 
tensities. The log-linear plots in FigJ5]reveal increasingly 
large deviations from the exponential behavior as we in- 
crease the pump intensity. All the distributions are well 
fitted by the stretched exponential function Ne~^ ci+C21 , 
1/ci providing a measure of the deviation from the ex- 
ponential function. The black line is the fitting function 
of the PDF with I p /Ith = 6.4. The fit is performed by 
setting the mean intensity and the variance of the fitting 
function equal to the experimental values. 

Note that an exponential statistics for the intensity 
corresponds to a Gaussian statistics for the field ampli- 
tude. Therefore, an exponential intensity PDF is charac- 
teristic of a speckles pattern, where each point receive the 
uncorrelated contributions of many uncoupled modes. At 
low pump, this is indeed the behavior displayed by the 
cavity field. However, when the pump increases, the in- 
creasing nonlinear coupling leads to a complex space-time 
dynamics and extreme events populate the tails of the 
PDF, providing a large deviation from Gaussianity. The 
L-shape of the statistics is clearly visible in the inset of 
Figl^l where the PDF are plotted in linear scale. 

To perform numerical simulations, we have considered 
the full model equations developed in Ref. 
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where no and n\ are, respectively, the amplitude of the 
homogeneous refractive index and the amplitude of the 
refractive index grating at the spatial frequency k c — k p , 
k c and k p being the optical wave numbers of the pump 
and cavity field, a is the nonlinear coefficient of the 
LCLV, and we have neglected the diffusion length due 
to elastic coupling in the liquid crystal. The dynamics 
of the liquid crystals is much slower than the settling of 
the cavity field, thus E c follows adiabatically the evolu- 
tion of hq and m. Then, by taking the wave propagation 
equation with the cavity boundary conditions, we obtain 
that 



E r = i 



[Ce ina J (2\n 1 \)} h Ce in °p-J 1 (2\n 1 \)E p , (2) 



where J m is the Bessel function of the first kind 
and of order to, and C is an operator accounting 
for the geometry of the cavity and losses, C — 
r 1 / 2 e i ' 5 S' :E e l(LoV i /2fe p ) e" 4(fc p' : 'i/ 2/) e l(LlV i /2fe p ) , with 1-T 
the photon losses, Lq = L — L\ and S the phase retar- 
dation in a round-trip. S x is a symmetry operator that 
inverts the x axis, thus accounting for the odd number 
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of mirrors, Vj_ is the transverse Laplacian and f± the 
position in the transverse (x, y) plane. 

For the numerical convenience, C has been imple- 
mented by using an equivalent nearly plane cavity with 
the y-axis inversion S y and two lenses just before and 



after the LCLV [161 ] ■ In order to account for the finite 



size of the LCLV and the diaphragm inserted near the 
lens, we have used in the simulations a diaphragm both 
in the space and in the Fourier plane. In particular, we 
have used a spatial filter with radius equal to 0.4 cm and 
a Fourier filter with radius k = 3.2 • 10~ 3 k p . For the 
other parameters, we have set the liquid crystal response 
time t = 0.1 s, the photon loss fraction 1 — T between 
0.7 and 0.8, and a is chosen in such a way that the in- 
tensity unit is the pump threshold for the activation of 
the cavity field. 
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FIG. 3: (color online) a) b) Numerical snapshots of the cavity 
field intensity at two different instant times; I p — 8, 1 — T = 
0.7; c) experimental snapshot of the cavity field, I p /Ith = 4.0. 
d) ID section of a 2D numerical intensity distribution. 

Starting with a random initial condition, we observe a 
transient speckle-like behavior, then a breaking of the 
y O — y inversion symmetry occurs and the coarse- 
grained intensity distribution \E c (x,y)\ 2 becomes non- 
homogeneous along the y axis. In Figj3] we report two 
snapshots of the cavity field intensity \E C \ 2 numerically 
calculated for I p = 8, 1 — T = 0.7 and displayed at dif- 
ferent instant times. In Fig[3^ we can note two large 
dark (bright) regions symmetrically distributed around 
the center, while in FigOb the left pair has disappeared 
and the other one is inverted. Typically, the cavity field 
shows a dynamical evolution with the succession of differ- 
ent spatial configurations, each living for a few seconds. 
For comparison, we display in Figj3t an experimental 
snapshot recorded for I p /Ith = 4.0. In Fig. 3d, we report 
a one-dimensional section of the numerical intensity dis- 
tribution at another instant time, exhibiting a very large 
narrow peak in a bright region. Its spatial size is larger 
than the experimental value because of a narrower spec- 
tral filtering that limits, for computational convenience, 
the number of numerical lattice points. 

In Figf5]we have plotted the numerical intensity \E C \ 2 
in the z — y plane. The yellow line at z = represents 
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FIG. 4: (color online) Cavity field intensity in the y — z plane; 
the yellow vertical line at z — indicates the position of the 
LCLV; I P = 10, 1 - T = 0.7. 



the position of the LCLV inside the cavity. In the plot, 
there is clear evidence of both the dynamical symmetry 
breaking of the intensity profile along y and of the y- 
axis inversion between the impinging and outgoing fields 
(at left and right of the yellow line, respectively) , consis- 
tent with the geometry of the experimental cavity, which 
is nearly spherical and made of an odd number of mir- 
rors. This inversion induces a nonlocal coupling between 
different spatial regions of the field, thus triggering the 
hypercycle amplification process that is essential for the 
symmetry breaking and the long tail statistics. Indeed, 
we have checked that simulations with a non-inverting 
nearly plane cavity generate an exponential statistics. 

As a consequence of the spatial symmetry breaking, 
the field exhibits large deviations from Gaussianity. The 
numerical PDF of the cavity field intensity are displayed 
in Fig J5] for different values of the pump. The inten- 
sities are rescaled in such a way that the PDF have 
the same value and slope at the origin. The distribu- 
tions are well fitted by the stretched exponential func- 
tion iVexp(— \Jc\ + C2I) (black dotted lines) and their 
tails are increasingly populated as the pump increases, 
in agreement with the experimental observations. The 
simulations show also that the symmetry breaking is 
stronger for high cavity losses. Moreover, for the used 
parameters, the refractive index uq is about four times 
less than ri\. The role of no is to adjust the global phase 
detuning S, but it has no effect on the PDF, as we have 
verified by removing it from the model. Thus, in the 
following we take n constant and 1 — T <C 1. 




FIG. 5: (color online) Numerical PDF of the cavity field in- 
tensity for 1 - T = 0.7 (a), 1 - T = 0.8 (b) and different 
values of the pump: I p /Ith = 6 (blue), 8 (red), 10 (green). 
The black straight and dotted lines are respectively exp(— I) 
and N exp( — \/ci + C2I), where N and d are set in order to 
fit the green distribution; (a) ci = 3.73, (b) ci = 2.83. 
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The y axis inversion (Fig|4| introduces a nonlocal cou- 
pling between different space domains with loops of am- 
plification, a mechanism similar to the hyper-cycle chain 
of reactions in the catalytic processes . Because of this 
mechanism, there is a sort of focusing, for which at some 
space locations the cavity field grows much more with 
respect to the surrounding places, giving rise to large 
amplitude peaks and, hence, to large deviations from the 
Gaussian statistics. The peaks last a few seconds, after 
that the hyper-cycle readjusts over new field configura- 
tions, yielding new peaks in other space locations. 

To elucidate this mechanism we derive a simple two- 
mode model, where only the evolution of the average re- 
fractive index m is kept, the average being performed 
over the transverse plane. For a nearly plane cav- 
ity, n\ satisfies approximatively the equation t^^- — 
—fix + IpFifii), where the first and second terms at the 
r.h.s. account respectively for the liquid crystal relax- 
ation and the grating feeding provided by the pump and 
cavity fields. As an ansatz for the function F , we take 
a cubic function, which describes a linear growth fol- 
lowed by saturation, due to multiple scattering and pump 
depletion [14J. 
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FIG. 6: The refractive index n+ as a function of I p . 

When the nearly plane cavity is replaced by a nearly 
spherical one, the above mean-field picture is modified by 
the nonlocal coupling due to the y o — y inversion. We 
thus have to consider two mean fields n + and n_, where 
the averages are performed over the upper and lower half- 
planes, respectively. Accounting for the inversion y — > 
—y of the cavity field, we have that, after a round trip, the 
grating h± is fed by the grating n=p at the opposite side, 
i.e, the equation for n is replaced by the two following 
ones for n±, 

r^ = -n ± + I p F(n T ), (3) 

where we have assumed implicitly that the cavity losses 
are very high, i.e., the cavity field is negligible after more 
than one round trip. 

The above equations are similar to a hypercycle model 
for two autocatalytic systems. The main features of F 



can be captured by the cubic function F(n±) — n± — n±. 
It is easy to show that for 1 < I p < 2, Eqs.Q have 
only one stable solution, n_ = n + = 0, whereas for 
I p > 2 a bifurcation occurs with the birth of the two 
stable asymmetric states n + = 2~3[1 ± (1^ — 4)5 /I p ]%, 

ri_ = 2~^[1 =F (I p - 4)i//j,]i, which break the y O -y 
symmetry, thus qualitatively accounting for the experi- 
mental and numerical observations. The bifurcation di- 
agram of n + is plotted in Fig|B] as function of the pump 
intensity I p . The experimental and numerical bifurcation 
intensity is about 4, rather than 2. This discrepancy is 
mainly due to the approximation of F(n±) with a cubic 
function. 

In conclusion, we have shown experimentally and nu- 
merically that extreme waves arise in a spatially extended 
nonlinear optical system. We have derived a simpli- 
fied mean-field model, showing that a spatial symmetry 
breaking is induced by a nonlocal coupling of the cavity 
field, which is responsible for a hyper-cycle type amplifi- 
cation of the field and leads to a non-Gaussian statistics. 
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